Magnetoconductance switching in an array of oval quantum dots 
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Employing oval shaped quantum billiards connected by quantum wires as the building blocks of 
a linear quantum dot array, we calculate the ballistic magnetoconductance in the linear response 
regime. Optimizing the geometry of the billiards, we aim at a maximal finite- over zero- field ratio 
of the magnetoconductance. This switching effect arises from a relative phase change of scattering 
states in the oval quantum dot through the applied magnetic field, which lifts a suppression of the 
transmission characteristic for a certain range of geometry parameters. It is shown that a sustainable 
switching ratio is reached for a very low field strength, which is multiplied by connecting only a 
second dot to the single one. The impact of disorder is addressed in the form of remote impurity 
scattering, which poses a temperature dependent lower bound for the switching ratio, showing that 
this effect should be readily observable in experiments. 

PACS numbers: 73.23.-b, 73.23.Ad, 75.47.-m 
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I. INTRODUCTION 

The ability to reduce the size of electronic circuits to 
the nanometric scale has lead to increasing interest in 
the properties of electron transport in the mesoscopic 
regime, and its dependence on externally tuned param- 
eters. Formation of two-dimensional (2D) structures 
of controllable geometry at semiconductor interfaces, so 
called electron billiards, sets the experimental grounds 
for investigating phase coherent transport of electrons. 
Open semiconductor quantum billiards serve as artifi- 
cial scatterers of highly tunable characteristics and have 
pioneered the understanding of the underlying physics, 
on both experimental and theoretical grounds. They 
are used to demonstrate and investigate a series of in- 
teresting phenomena on the mesoscopic level, such as 
shot noise in transport through charged dot^*^, Fano 
resonances^i^i^, Andreev tunneling and reffection^i^i^, 
decoherence in ballistic nanostructures^^^^^ as well as 
classical to quantum transitions and imprints of nonlinear 
dynamics^ii^ii^ii^. Further, the geometry of a conducting 
structure is shown to have a major impact on the result- 
ing transport phenomena^^^^^. The magnetoconductance 
of such nanodevices proves as an essential signature for 
the underlying interference phenomena and has there- 
fore been studied extensivel}siiii^ii^i2P . The Aharonov- 
Bohm effectp^ is directly observed in systems of quan- 
tum ringsi^i2^i2^i2i, but also plays a central role in de- 
scribing magnetoconductance fluctuations in more com- 
plex mesoscopic systems in weak magnetic fields^^^^^. At 
higher magnetic field strengths the quantum Hall effect 
sets in, accounting for a steplike varying magnetoconduc- 
tance, formation of edge states and characteristic multi- 
channel fluctuations in the transmission spectra3ii2^i2^. 
Localization effects and conductance fluctuations mani- 
fest themselves in a large variety of open quantum dot 
systems, regardless of whether ballisti o-'^^i-'^^i^^i^^i^-'^i^^ or 
diffusive —i^i^i^ transport is considered. Assembling in- 



dividual dots into coupled arrays or lattices gives rise to 
new features of the system's overall response, depend- 
ing on the type and strength of couplin g^^i^^i^^i^^i^^i^^ . 
Of particular interest are systems where the interplay 
between the various effects of electron transport men- 
tioned above can be used to achieve a tunable quan- 
tum conductance, in terms of designing the size, shape 
and material specific features of the conducting device, 
as well as varying macroscopically accessible parame- 
ters such as externally applied fields, temperature, and 

gate voltages controlling the coupling strength between 
constituent^a^^^^^^2 

In this article we exploit the dependence of the con- 
ductance on the geometry of a 2D electron billiard and 
examine its functionality as a switch when a magnetic 
field is turned on. Employing oval shaped billiards as 
the building blocks of a linear quantum dot array, we 
aim at a maximal finite- over zero-field ratio of the con- 
ductance by optimizing the system within an achievable 
parameter range. The switching effect arises from the 
phase changing effect of the applied field, which raises a 
suppression of transmission present for a certain deforma- 
tion of the oval. The assembly of dots into a chain even- 
tually leads to banded transmission spectra for a large 
number of dots, with details depending on the interdot 
lead length. The conductance, taken as the thermally 
averaged transmission function, oscillates with increas- 
ing field strength; at higher fields edge states form, which 
conduct ideally. The switching ratio corresponding to the 
first magnetoconductance maximum acquires a multiple 
value by adding one more oval to the single one, while it 
fluctuates for further added dots. The impact of impuri- 
ties may enhance or weaken the switching effect, whether 
or not they block the leads coupled to the dots, imposing 
a temperature dependent lower bound on the switching 
ratio in the presence of weak disorder. 

In section II the setup and geometry of the 2D struc- 
ture are specified and the theoretical framework as well 
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as computational approach are presented. In section III 
the main features of the obtained transmission spectra 
are discussed, along with a description of the underly- 
ing mechanisms. This is followed by an analysis of the 
switching ratio in dependence of the deformation of the 
billiard shape, the magnetic field strength and the length 
of the multidot chain at different temperatures, in order 
to determine a device setup optimal for switching, within 
an achievable parameter range. Finally, the modification 
of the switching ratio in the presence of disorder is stud- 
ied. Section IV provides a summary of results, concluding 
on the functionality of the switching mechanism. 



II. SETUP AND COMPUTATIONAL 
APPROACH 

The confining potential of the single dot is assumed to 
be of hard wall character, leading to Dirichlet boundary 
conditions for the wave function. We use an oval billiard, 
whose shape is parametrized as^ 
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where (j) G [0,27r]. The parameter S tunes the defor- 
mation of the dot, which becomes a circular billiard of 
radius R if S = 0. In this case the classical dynamics of 
the closed system is integrable, whereas for 6 > it be- 
comes non-integrable with mixed phase space^^iM. For 
reference with respect to the device specific parameters, 
a mesoscopic size of R = 220 nm is employed. At the 
right and left ends of the elongated structure semi-infinite 
leads of width W = 0.3 R are connected, representing the 
coupling to electron reservoirs. The use of semi-infinite 
leads models the ideal case of vanishing refiection of the 
electrons upon reaching the reservoirs. In the multidot 
case this single cavity is replaced by a chain of N iden- 
tical oval dots connected to each other through leads of 
length L, where L equals the distance between adjacent 
oval edges, aligned with the semi-infinite leads possess- 
ing the same width. Fig. [1] provides a picture of the 2D 
structure for A/" = 2 connected dots. 

Restricting ourselves to low temperatures and a small 
system size we neglect inelastic processes, and do not 
account for electron-electron or electron-phonon interac- 
tions. The single particle Hamiltonian is, within an ef- 
fective mass approach, of the form 



H 



2meff 



V{r) 



(2) 



where we choose a value of meff = 0.069 rrie corre- 
sponding to GaAs, rrie denoting the electron mass. V{t) 
is the hard wall potential, and the vector potential A 
produces a magnetic field perpendicular to the plane of 
the structure, over which it is homogeneously extended 
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FIG. 1: Geometry of the billiard forAr = 2,5 = 0.5,L = 
W — 0.3 i?, consisting of the repeated module and the outer 
leads which represent the connection to reservoirs (see text). 



with strength 5, dropping off linearly to zero in the ex- 
terior leads. We will concentrate on the magnetocon- 
ductance switching effect at a very low magnetic field 
strength (~ 0.02 T), where the Zeeman splitting for 
GaAs (~ 3.6 /ieV) is negligible (~ 0.1%) with respect 
to the Fermi energies we consider, and therefore do not 
take into account the coupling of the electronic spin to 
the magnetic field. The Hamiltonian is discretized on a 
tight-binding lattice, with the magnetic vector potential 
incorporated through Peierls' substitution. The coupling 
of the system to the external semi-infinite leads placed on 
the left (/) and right (r) of the billiard is described by self- 
energies which are analytically obtained for 5 = 
and contribute non-Hermitian blocks to the Hamiltonian 
matrix. From the single-particle Green's function of the 
system 



G(E) = (H + E^ + E^ 



(3) 



the part G^^ describing the propagation from the left to 
the right lead is computed using a parallel implementa- 
tion of the recursive Green's function method (RGM), 
where a decomposition scheme among communicating 
processors allow for the computation to be done in a par- 
allel manner^^ . In the multidot case the chain is built up 
by a repeated module, which consists of the oval cavity 
with leadstubs of length L/2 on the right and left (see 
Fig. [1]). Having found G^z for one module, we calculate 
the Green's function connecting the two outer leads using 
a modular variant of the RGM, which was originally pre- 
sented in Ref. 0. In this algorithm the Green's function 
of the joined module is calculated using the Dyson equa- 
tion. The transmission of the device is finally evaluated 
via the Fisher-Lee relations"^^, T{E) = Tr (P^GP/Gt), 
with Vijj. = — ^l/r)' worthwhile noting that 

in the two-terminal device we encounter, even in the pres- 
ence of a magnetic field, the transmission function is sym- 
metric under the exchange of the contact leads, i.e. the 
transmission from left to right equals that from right to 
leflj^. The computed Green's function of the system is 
also used to calculate the local density of states (DoS) 
at site r through the relation p{r,E) = (r|A(£^)|r)/27r, 
where A = GPG^ is the spectral function and P gener- 
ally a weighted sum of P^ and P^ according to the Fermi 
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distributions of incoming states in the two leads. In the 
cases presented here, we have chosen F = F^, i.e. p(r, E) 
corresponds to the probabihty density resulting from an 
incoming monochromatic wave of energy E from the left 
lead. 

The calculated transmission determines the macro- 
scopically measurable conductance of the device. In the 
linear response regime at low temperature the conduc- 
tance for given Fermi energy Ep can be obtained by the 
Landauer formulai^-'— : 



h J-oo 



T{E)Fth{EF;E)dE 



(4) 



with 



Fth{EF;E) 



df{EF;E) 
dE 

sech? 



1 



E-Ef 
2kBe 



(5) 



where f^Ep^E) is the Fermi distribution function cen- 
tred around Ep^ and thus it essentially equals the ther- 
mally averaged transmission around the electron Fermi 
energy, with a width determined by the temperature O. 



III. RESULTS 

A. Transmission spectra 

Concentrating on the transmission in the deep quan- 
tum regime, we restrict the energy of the incoming elec- 
trons such that only the transversal ground state of the 
leads is energetically available. Thus the dimensionless 
channel number k = kW/n^ where k = \/2mei^E/h^ takes 
on values in the range 1 < < 2, that is, in the first 
channel of transmission. For the size of the device spec- 
ified, this corresponds to a Fermi energy in the range 
1.2 meV < Ep < 5 meV. A detailed analysis of the trans- 
mission within the first channel in terms of the quantum 
states in the single oval dot has been presented in Ref. 
Here we focus on the modification of the transmission 
when dots are connected to form an array, as well as 
the conductance of the device (where the details of the 
transmission are thermally averaged out) as a function of 
the geometry parameters and the magnetic field strength. 
The zero- and finite-field transmission T(/^) through the 
device is shown in Fig. [2] for different numbers of dots 
N in the chain, with deformation parameter S = 0.5 and 
inter dot distance L = W. As the channel number Hi 
measures the wave number in units of tt/W, T{hz) de- 
pends only on the ratio W/R. Our calculations show 
that changing W/R within 0.2 < W/R < 0.4 introduces 
mainly a shift in T(/^) according to the implicit energy 
scaling, i.e. the transmission is largely determined by the 
geometry of the billiard and not by the leadwidth. For 
values of W/R > 2, the transmission obviously has to ac- 
quire the value of the unperturbed quantum wire. In the 



following we restrict ourselves to the case of W/R = 0.3. 
The zero-field transmission in the single dot case consists 
of a rather smoothly varying background, on which sharp 
Fano resonances are superimposed. In the multi-oval case 
these sharp resonances are TV-fold split (this very small 
splitting is generally not resolved on the scale of Fig. [2]), 
while, additionally, Breit-Wigner (BW) type resonances 
of varying width emerge, firstly for = 2, and subse- 
quently undergo a splitting into A — 1 sub-peaks for an 
array of A dots. For sufficiently many dots (represented 
in Fig. [2] by the case of A = 20), the multiply split res- 
onances saturate into bands of densely positioned peaks, 
which is reminiscent of the band structure of energy lev- 
els in a periodic quantum system. In the presence of the 
weak field the smooth background transmission is overall 
increased, the sharp resonances are slightly shifted in en- 
ergy and the transmittive bands for large A are broader. 

The sharp Fano resonances of the single dot case cor- 
respond to quasibound states that are strongly localized 
within the dot, at energies which coincide with eigenen- 
ergies of the closed oval billiard without attached wires. 
There is a series of equidistant resonances correspond- 
ing to the different longitudinal modes for each given ex- 
cited transversal mode in the oval^^. The spatial confine- 
ment of such states decouples them from the leads and 
consequently they do not contribute significantly to the 
conductance of the device. The states that can provide 
a substantial contribution to the conductance are those 
with a longitudinal spatial extension unto the openings of 
the leads. They are strongly coupled to the leads of the 
open system in the case of constructive interference at 
the openings, leading to a broad transmission maximum. 
We refer to these states, which extend from the oval into 
the leads, as leaking states. The number of leaking states 
is determined by the allowed excitations inside the cavity 
subject to the constraint of the energy being within the 
first channel. The interference of leaking states belong- 
ing to different transversal modes generates broad humps 
in the single dot transmission, where the transmission is 
substantial over a finite energy interval (constructive in- 
terference), separated by points of vanishing transmission 
(destructive interference). The slowly varying envelope 
behavior of the transmission spectrum exhibits a wide 
energy range where the overall transmission is strongly 
suppressed. For the specific shape of the cavity corre- 
sponding to the chosen value of S = 0.5, this suppression 
valley is centered around the middle of the first channel. 

In order to analyze the transmission of the multidot 
chain, in Fig. [3] (A) we focus on the transmission around 
the BW resonance appearing for N = 2 dit k = ^ 
1.384, and show its (A — l)-fold splitting for increas- 
ing A. Also the sharp Fano resonance just below is in- 
cluded, whose splitting (of the order of Ahc ~ 2 x 10~^ 
or AE ~ 0.1 fieV) remains unresolved even at this scale. 
The A-fold splitting of the Fano resonances is a con- 
sequence of the degeneracy of the confined single dot 
eigenstates in the case of A dots, which are coupled 
very weakly through the connecting lead due to their 
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FIG. 2: (Color online) Transmission spectra in the first transversal channel for varying number of dots N with deformation 
parameter (5 = 0.5 and connecting lead length L = VK = 0.3i?, at5 = (solid black line) and B — Be ^ 2^ mT (dotted red 
line) . 



strong localization within the ovals. It is thus similar to 
the splitting of the energy levels of atoms brought to- 
gether to form a weakly bound molecule, with an energy 
split proportional to the interatomic coupling^^. The 
BW type resonances of the multidot case, which are nar- 
rower (wider) at energies where the single dot transmis- 
sion T^^^^\k) is lower (higher), are of different origin: 
They arise from the resonant tunneling of the incoming 
wave through the system of the ovals and the connect- 
ing bridges. Indeed, the emergence of these resonances 
and their (A^ — l)-fold splitting can effectively be deduced 
from the ID scattering through potential barriers (or 
equivalently, N — 1 resonators), where the transmission 
amplitude of scattering through each barrier possesses 
an energy dependent norm and phase. Two barriers a, 
13 with transmissions Tq,, T/3 give the total transmission 



fa/3 



(6) 



where Ra//3 = 1 — ^q;//3 and is the phase shift ac- 
quired by reflection from P to a and back to p. For 
Ta = Tp = T^^^^^ and cc k, this gives rise to resonance 
peaks in T^^^^^ which are equidistant in k, and have a 
width that increases with T'^^\ In our case though, due 
to the structure of the ovals that constitute the barriers, 
the phase shift is not linear in Hi. This perturbs the pe- 
riodicity of the resonances, as we observe for T^'^\ti) in 
Fig. [21 or equivalently, yields an energy dependent effec- 
tive resonator length L{i<i) oc 0{i<i)/ n. Formula (|6]) can be 
iterated to obtain the transmission ior N ^ 2 ovals, i.e. 



TW = T^p{ T« = TW,T0 = r(^-i);e«,;3 = ^1,^-1 ), 
where 9 now results from reflections between 1 and N —1 
barriers. The (7V-l)-fold splitting of the T^^) resonance, 
shown in Fig. [3] (A), and the saturation into a band in 
the transmission spectrum for large N ^ are then repro- 
duced for a system that is symmetric under the exchange 
a ^ P (which, in our case, renders the dots identical), 
provided that the phase difference between transmission 
and reflection amplitude of the single barrier is equal to 
±7r/2, as is the case for the single oval with symmetric 
leads. Varying the resonator length modifies the condi- 
tions for resonant transmission by shifting the resonances 
in energy and changing their periodicity. In Fig. [3] (B) 
the transmission through N = 2 connected dots, as well 
as the normalized conductance at = 0.2 K, are plot- 
ted over the energy range of a single dot transmission 
hump, for varying connecting bridge length L. With a 
slight increase in L {L/W = 1.0,1.1,1.2,1.3) the BW 
resonances are shifted to lower energy, and for longer 
bridges {L/W = 20, 100) the number of resonances in the 
same interval increases. We notice that the center posi- 
tions of the (split) Fano resonances are unaffected by the 
variation of the bridge length. Detailed features of the 
transmission lineshape, such as the Fano resonances and 
the BW resonance peaks for large L, are washed out by 
thermal averaging, making their contribution to the con- 
ductance negligible compared to the smooth background. 

As we see in Fig. [3] (A), the addition of a dot to the ex- 
isting chain at resonance energy, lowers the transmission 
from unity to the single oval value T'^^^hc) at that energy 
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FIG. 3: (Color online) (A): T{k) for B = (solid black line) 
and B — Be (dotted red line) for L — W and varying N ^ 
in the vicinity of Kp ^ 1.384, with labels a, b, c, d, e for 
the resonances referred to in the text, (B): T{Hi) (solid black 
line) and g{K, ; B = 0.2 K) (dashed magenta Une) for a single 
dot (bottom) and for two dots with varying bridge length L, 
within a small window of the channel number k covering the 
energy range of a single dot smooth hump. 



(note that the transmission at the dips between the res- 
onances can acquire values even lower than T'^^^i^) ). 
In particular, the transmission at the energy position 
of the central resonance at />: = oscillates between 
unity and T^^\Kp) with even and odd respectively: 
T^^''''''''\kp) = 1, T(^^^^)(/^p) = T^^\np). Furthermore, 
the resonances for each N are positioned symmetrically 
around Kjp , SO that the forming bands in the transmission 
for large N are centered around the T^'^\k,) resonance 
peaks. In Fig. 2] this behavior of the transmission func- 
tion for varying number of dots is illustrated in terms 
of the states forming in the system for TV = 1,2,3,4 
dots, by plotting the zero-field local DoS at the energies 
(rows a,b,c,d,e) of the resonance peaks labeled (with cor- 
responding letters) in Fig. [3] (A), for an electron incident 
on the left. The spatial oscillations of the DoS in the 
incoming lead come from the interference of the incom- 
ing wave with the wave that is backscattered from the 
dot array. Their absence is a signature of a resonance 
peak in the transmission spectrum, as there is no over- 
all backscattering and transmission is unity. It must be 
noted that the colormap for the DoS in each of the sub- 
plots is normalized to the maximal value, such that same 
colors at different sub-plots do not represent equal abso- 
lute values (which are irrelevant in the present analysis). 
Starting with the single oval in the first column [Fig. [H 



(al)-(el)], we see that the incoming wave is reflected at 
all energies (a) to (e) around A^p, leading to a transmission 
significantly less than unity {T^^\K,p) ^ 0.07). When a 
second oval is added, there is an energy between (a) and 
(e), namely n = represented by row (c), for which the 
backscattering of the single oval is cancelled by the pres- 
ence of the added oval and the connecting bridge: the 
wave is multiply reflected between the two ovals through 
the bridge, resulting in a quasi-standing wave along the 
chain [see (c2)] that constitutes a resonant state for the 
open system, leading to a transmission of unity. A third 
oval added in front of the two introduces the backscatter- 
ing again at np [see (c3)], while the transmitted part from 
this first oval is perfectly propagated through the remain- 
ing two as in (c2). Thus the transmission in (c3) equals 
the single dot transmission, T^^\tip) = T^^^^Kp). The 
backscattering of the third oval at Kp is cancelled by ad- 
dition of a fourth oval [see (c4)], just as we went from (cl) 
to (c2), so that the T^'^\i<ip) resonance peak is recovered 
in T^^\i<ip)^ although with a smaller width. Thus, the se- 
quential cancellation of the backscattered wave leads to 
the even-odd oscillations of T^^\hip) seen in Fig. [3] (A). 
Resonant states are also accessed for 3 and 4 dots in (b3) 
and (a4) below hZp and in (d3) and (e4) symmetrically 
above tip. Similarly, for each number of dots N there 
are N — 1 accessible resonant states, including the one at 
tvp for even TV, at energies symmetrically positioned with 
respect to K,p. Just as the T^^^ resonance is recovered 
m tW, each T(^) resonance is recovered at multiples 
of where the resonant state in the chain can be de- 
composed into multiple connected resonant states. We 
notice that the two branches of resonances, one below 
and one above hZp^ are associated with two different leak- 
ing eigenstates of the single oval with closed leads - they 
inhabit, for example, the central oval in (b3) and (d3), 
respectively. Their interference in the open single oval 
system forms the scattering wave in column (1). These 
three wavepatterns are combined among the N ovals in 
the open chain, to form the — 1 resonant states lead- 
ing to the peaks around f<ip. The formation of resonant 
states occurs similarly around all T^^^ resonances of BW 
type (seen in Fig. [2]). Characteristically, moving from 
a T^^^ resonance to the next one at higher energy adds 
a node in the quasi-standing wave within the two ovals 
and the connecting bridge. Increasing the length of the 
bridge shifts the resonances to lower energies and reduces 
the /^-distance between them, as the wavelength in the 
quasi- standing wave overall increases, in accordance to 
the effective resonator picture described above. 

Conclusively, there are two types of resonances to be 
distinguished in the transmission spectra for the array 
of N dots: {i) the series of equidistant Fano resonances, 
arising from the confined single dot excitation modes in 
the continuum of the channel, which are A^-fold split 
due to coupling between the ovals, and (n) the series 
of nonequidistant BW resonances, resulting from reso- 
nant tunneling states that form in the chain, which are 
{N - l)-fold split. 
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FIG. 4: (Color online) Zero- field local DoS for N — 1, 2, 3, 4 dots with 6 — 0.5 and L — W for energies in the vicinity of the 
(N — l)-fold split resonance peak ai k = K,p ^ 1.384, with incoming electron on the left. Rows (a), (b), (c), (d), (e) correspond 
to the energies of the resonances labeled with the same letters in Fig. [3] (A). The colormap for the DoS in each sub-plot is 
normalized to its maximal value, further the color maps the square root of the DoS to enhance contrast. 



Following the discussion above, we now consider the 
impact of the perpendicular homogeneous magnetic field 
on the transport through the device. When the field is 
switched on the phases of the different states forming in 
the ovals are modulated, and consequently the interfer- 
ence of the states contributing to transmission changes. 
Thus, depending on the field strength, the transmission 
spectra for the single and multiple dots are accordingly 
modified. As we see in Fig. O the weak field of 20 mT 
introduces dramatic changes in the spectra. The slowly 
varying background of the single oval case is generally 
raised throughout the channel, removing the character- 
istic suppression around its middle in the absence of the 
field. The overall very high transmission is interrupted 
by series of dips in its lineshape. The sharp Fano res- 
onances undergo only a very slight energy shift (visible 
for the Fano resonance in Fig. [3] (A) ), because the spa- 
tial distribution of the wave function remains practically 
unaffected by the low field chosen. Again the multidot 
chain provides a more complex transmission spectrum, 
resulting from the subsequent matching conditions for 
the wave function at the connections between the dots. 
The BW and Fano resonances are multiply split like in 
the field free case, and dips and plateaus become sharper 
and more pronounced as dots are added to the chain, 
saturating into a banded transmission. In contrast to 
the field free case, the transmission pattern is now dom- 
inated by narrower gaps and wider transmittive bands. 
Thus, also for the long chain of dots the overall trans- 
mission is drastically raised by the applied field. A more 
detailed analysis of the modification of the conductance 
with varying field will be presented in the next subsec- 



tion. 



B. Conductance switching 

The normalized conductance g = of the chain of 

oval dots is shown in Fig.[3](B) (dashed curve) for a tem- 
perature G = 0.2 K over a part of the first transmission 
channel, in direct comparison to the transmission func- 
tion. The parameter hi now represents the scaled Fermi 
energy of the incoming electrons, around which the trans- 
mission function is thermally averaged. At zero temper- 
ature conductance and transmission are equal, but as O 
is increased peaks and dips in the spectrum become less 
pronounced due to the increased range of contributing 
energies. As mentioned above, already at the low temper- 
ature chosen, the detailed structure of the transmission 
is essentially lost: the sharp resonant peaks are washed 
out, reflecting their negligible contribution to the conduc- 
tance. Also the formation of sharp transmittive bands 
for the multidot chain is relaxed with thermal averaging. 
For long interdot leads (Fig. [3]^d),(e)) the conductance 
features follow the trend of the single dot case, that is, it 
exhibits similar humps in energy, yet with smaller ampli- 
tude. Similar modifications of the transmission spectra 
through thermal averaging hold for the conductance pro- 
file in the presence of the magnetic field. 

A key feature of the oval shaped cavity is the formation 
of the wide suppression valley in the transmission spec- 
trum of the first transversal channel, which is essentially 
retained also for the conductance at low temperature. In 
order to demonstrate the suitability of the chain of dots 



as a magnetically induced conductance switch, we exploit 
the lifting of this suppression when the field is turned on, 
aiming at a high ratio of finite- over zero-field conduc- 
tance. In the following we optimize the switching ratio 
taking into account all relevant parameters (^, 5,L,A^), 
as well as finite temperature and impurity scattering ef- 
fects (see subsection III.C). First we consider the quan- 
tity G™ which is the zero-field finite temperature con- 
ductance minimized with respect to the position of the 
Fermi energy in the first channel. 

In Fig. El = 2t^GX is plotted as a function of 
S at different temperatures for a single oval dot. We see 
that an optimal value for is obtained around S = 0.5, 
with a small dip at ^ = 0.55, while it increases for larger 
or smaller deformation of the oval. It must be noted here 
that the modification in the spatial extension of the oval 
for a change AS ^ 0.05 is of the order of 1%, a chal- 
lenging accuracy for an experimental realization of the 
device. We therefore keep the roughly optimized value 
of J = 0.5 as a reference for the following analysis. As 
shown in the inset of Fig. [5l the channel number k,^^^ 
of this minimum depends approximately linearly on 6 
where the corresponding Fermi energies are located close 
to the center of the first channel. For S = 0.5 we have 
i^c = /^™(^ = 0.5) ^ 1.46 in the single dot case (TV = 1). 
This shift of the optimal Fermi energy, that holds for all 
temperatures considered, is due to the modification of 
transversal modes inside the dot, which are shifted to 
higher energies as the oval becomes narrower with in- 
creasing S. 

The single dot switching ratio S^^^^\B) = 

Gin~^\B) / G^^~^^ dit = tZc is shown in Fig. [6] for vary- 
ing magnetic field strength at different temperatures. As 
S^^\B) equals the finite field conductance normalized 




0.35 0.4 0.45 0.5 0.55 0.6 0.65 
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FIG. 5: Minimal zero-field conductance (see text) as 

a function of the deformation parameter S for a single dot, 
at temperatures (bottom to top) B = 1.0, 1.1, 2.0 K; the 
inset shows the change of the optimized channel number k™^^ 
with (5 at G = 2 K (the dependence is the same for the other 
temperature values). 
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FIG. 6: Single dot switching ratio S^'^^ at n = Kc ^ 1-46 
as a function of the magnetic field, for temperatures (top to 
bottom) B = 0.7, 1.0, 1.4 K; the inset shows the irregular 
oscillations for low field strengths. 



to G^fj, it describes the changes of the conductance in- 
duced by the field. For low field strengths (inset of Fig.[6|) 
the modulation of the phase of the longitudinal states in 
the dot leads to Aharonov-Bohm (AB) like oscillations 
in the conductance. At the energies we consider here, 
only three of these leaking states are present How- 
ever, the presence of more than two channels inside the 
dot gives rise to the superposition of magnetoconduc- 
tance oscillations, so that S'^^^B) loses the periodicity 
expected for AB oscillations of a ID quantum ring. As 
the field strength is increased, apart from their phase, 
also the spatial distribution of the states in the dot is 
affected. Confined states are eventually deformed into 
leaking ones, opening further channels for the transmis- 
sion. The first magnetoconductance peak at Be ^ 0.02 T 
is seen to be the highest in the low field regime, giving a 
switching ratio of S^^\B = Be) 65 at 8 = 0.7 K. For 
higher field strengths the transmittive states are gradu- 
ally localized into edge states (with Larmor radius ^ R/4 
for 5 > 0.8 T) along the border of the cavity, all within 
the first magnetic Landau level^^. Following the edges 
of the billiard, the electrons are now more easily trans- 
mitted, resulting in an increased overall conductance. At 
a field strength of 5 ~ 1.2 T these modes become per- 
fectly transmittive along the edges of the structure, and 
the switching ratio reaches a plateau of maximal value. 
For even higher magnetic field strength the transmission 
decreases drastically as the incoming electrons gradually 
fail to overcome the magnetic barrier provided by the 
first Landau level, and the conductance drops to zero. 
At higher temperatures the features of the magnetocon- 
ductance remain; however, as a broader energy window 
with higher transmission parts is contributing to the ther- 
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mal averaging, the switching ratio is generahy lowered, 
because G^J^ increases. Also the amplitude of the oscil- 
lations decreases with temperature, as the magnetically 
induced changes in the detailed structure of the transmis- 
sion have a smaller impact on average. For N > 1 the 
magnetoconductance behaves similarly, but the switch- 
ing ratio overall acquires higher values, because of the 
even lower zero-field conductance, resulting from the for- 
mation of gaps in the transmission spectra. 

The magnetoconductance is calculated for spinless par- 
ticles and hence does not describe electronic transport for 
high magnetic field strengths. But, as we are aiming at a 
high switching ratio, we concentrate in the following on 
the first maximum occurs 
approximately at the same field strength Be ^ 20 mT for 
all considered numbers of dots N. For this weak magnetic 
field we can neglect the Zeeman splitting. In Fig. [71 Sc^^ 
is presented for a varying number of dots in the chain, 
again at different temperatures. We allow for the param- 
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FIG. 7: Switching ratio Si^^ for varying number of dots N 
with connecting bridge length L — W ^ for different tempera- 
tures G. 

eter ni^^ , which represents the scaled Fermi energy of the 
incoming electrons, to be optimized individually to min- 
imize G^^^ for each number N . At sufficiently low tem- 
perature, by connecting a second oval to the single one we 
gain a substantial factor with respect to the increase from 
S^c^ to Sf^ ( ^ 320 for 6 = 0.7 K), which, as pointed 
out, results from the lower zero-field conductance. For 
N > 2 the switching ratio fluctuates around a temper- 
ature dependent mean value, due to its high sensitivity 
with respect to the optimized G^^^ at low temperatures, 
which changes for each N. At higher temperatures the 
fluctuations are weakened, but Sc^^ is then also lowered 
drastically. 

It is obvious that the optimization of the switching ra- 
tio strongly depends on the temperature: High switching 
ratios require low temperatures, < 2 K for our setup. 
Nevertheless, we see that the current switching function- 
ality of the device is signiflcantly enhanced throughout 



the temperature range considered, by taking e.g. two 
dots instead of a single one. 



C. The impact of impurities 

Let us explore the impact of impurity scattering, i.e. 
disorder, on the magnetoconductance. This is imple- 
mented in the form of remote impurity scattering in the 
presence of a modulation-doped layer above the 2D struc- 
ture. We consider pointlike negatively charged impurities 
of 2D density nimp distributed on a plane at distance d 
above the 2D electron gas (2DEG), excluding them from 
the region of the semi-inflnite leads. The plane is par- 
titioned into small pieces of area l/nimp, within each of 
which one impurity is placed at random position, thus 
constituting a quasi-random distribution of impurities, 
with an upper bound on their local concentration. The 
electrostatic potential of each impurity is screened by the 
2DEG at the plane of the device structure, so that the ef- 
fective potential that an electron feels at distance r from 
the impurity is modeled by^^ 



where 



A{d) = 



<7TF(l + gTF^) 

47reoe5 



^TF 



(7) 



(8) 



with 65 denoting the relative permittivity of the material. 
The Thomas-Fermi screening wave number (^tf is, for 
the low temperatures considered, approximated by (^tf ^ 
2/aB, where is the effective Bohr radius. As typical 
values for a GaAs semiconductor we take = 13.8 and 
aB = 9.8 nm. 

As the distance d of the impurity layer is made very 
short {d < 30 nm in the present scaling), the correspond- 
ing transmission spectra (not shown here) are drastically 
changed with respect to the clean case (see Fig. [2j), as a 
result of the inffuence of the impurity potential on the 
transport through the device. The randomized potential 
landscape in the dot chain leads to a spatial deforma- 
tion of the existing states and a breaking of the symme- 
tries present in the clean system: The sharp Fano res- 
onances are shifted due to the perturbation of the con- 
fined eigenstates in each dot, differently for each individ- 
ual impurity configuration. The impurity potential also 
changes the energies of the leaking states, which results 
in modified conditions for their coupling to the leads, 
so that the broad transmission maxima are shifted, too. 
Additionally, new transmission peaks are introduced by 
leaking states that did not contribute in the clean case 
due to their symmetry^^. For not too short impurity 
layer distance though, the described suppression valley 
in the conductance of the clean system is retained, still 
making it sensible to speak about magnetic conductance 
switching. The effects of disorder are of course enhanced 
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FIG. 8: Configuration mean and standard deviation of the 
switching ratio S^^^ of two connected dots as a function of 
the impurity layer distance d for different temperatures. The 
dashed lines give the values of the disorder free case. 



with increasing impurity density; we use here a value of 
^imp = 0.0025 nm~^. This rather high density of remote 
impurities is employed here in order to intensify their im- 
pact on transport in our simulations, whereas in practice 
cleaner samples are realizable for use in semiconductor 
nanostructures^^^^^ . 

In Fig. [8] the switching ratio is shown as a function of 
the distance d from the impurity layer for two connected 

(2) 

ovals. The values of S^^^ for each d are the average over 
27 configurations of the randomly distributed impurities. 
When the impurity layer is closer to the 2D conducting 
structure, the average switching ratio is in general lower 
than its value in the clean system, the latter being prac- 
tically reached for a distance d > do{&), depending on 
the temperature. For O 1 K we have do ~ 100 nm, 
corresponding to a transport mean free path Itr ^ 24 fim 
in the first Born approximation^^. Nevertheless, the rel- 
atively large deviations from the mean indicate that, for 
each d < do{Q)^ there are certain impurity configurations 
that provide a switching ratio much higher or lower than 
the average. This is due to the high sensitivity of Goff 
with respect to the potential pattern that is formed on 
the plane of the array. If the impurity configuration is, for 
example, such that a potential maximum is blocking the 
opening of a cavity to a lead, then Goff is suppressed, as 
the wave coming from the lead is strongly backscattered. 
This backscattering can be lifted when the magnetic field 
is turned on, leading to an overall increased switching ra- 
tio for this configuration. On the other hand, when the 
configuration of the impurities does not block the leads, 
Goff in the suppression valley is slightly higher compared 
to the clean case due to the additional resonances in the 
transmission, causing a reduced switching ratio. Thus, 
at distances where the potential on the 2DEG plane is 
not too strong to permit transmission at all, the ran- 
domly distributed impurities lie within a broad varia- 
tion between the cases of blocking and non-blocking con- 



figurations, keeping the deviations from the mean high. 
When the impurities are put too close to the 2D struc- 
ture {d < 30 nm), the shape specific suppression feature 
of the zero-field transmission is essentially lost, so that 
the overall conductance is practically unaffected by the 
field strength, which thus minimizes the switching effect. 

(2) 

For larger impurity layer distances the mean S^^^ even- 
tually saturates into the clean case value with decreasing 
deviations, as the potential becomes too weak to affect 
the transmittive states in the dots. 

Using random impurity distributions to investigate the 
functionality of magnetic current switching in a more re- 
alistic environment, one can speak of a temperature de- 
pendent lower bound of the switching ratio (see Fig. [8|) 
depending on the specific setup. This lower bound is in- 
creased as the influence of disorder is suppressed, that 
is, when a longer mean free path for the electrons is 
achieved. Technological progress actually makes it feasi- 
ble to reach mean free paths in heterostructures compara- 
ble to the size of realizable nanoscale devices —i^*^. The 
almost ballistic nature of electron transport then allows 
for controllable conductance switching at low tempera- 
tures, in the sense that it is determined by the specific 
shape of the conducting device, the electron energy and 
the applied magnetic field. 



IV. CONCLUSIONS 

Having investigated the transmission properties of a 
linear array of equidistant identical oval shaped quan- 
tum dots, we demonstrated the functionality of such a 
structure as a magnetically controlled switching device 
in the deep quantum regime. The switching effect arises 
from the lifting of a deformation specific suppression in 
the transmission of the oval when a weak perpendicular 
field is turned on. The suppression valley in the trans- 
mission results from the destructive interference of states 
in the dots that are strongly coupled to the leads, and 
is specific to the elongated shape of the single billiard. 
This makes the effect relevant in systems of similarly 
shaped dots (e.g. elliptical). The switching ratio oscil- 
lates with the magnetic field strength, but as the effect 
is prominently present even at very weak fields, we have 
concentrated on its first peak. We have shown that the 
extension of the single dot into a chain of dots causes a 
much higher switching ratio, due to a stronger suppres- 
sion of the zero-field conductance. However, we point 
out that almost optimal switching can be obtained by 
connecting only one more dot to the single one, giving a 
multiple value for the switching ratio while keeping the 
system size small. This could make the device practi- 
cally advantageous but also favors quantum coherence 
itself, which is the principal requirement for the inter- 
ference effects to take place. The efficiency of switching 
is lowered with increasing temperature, as the desired 
shape specific characteristics of the transmission spectra 
are thermally washed out, which poses a limitation to 
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low temperatures (up to about 2 Kelvin). In spite of the 
possibility to achieve mean free paths of the 2DEG much 
longer than the extent of the studied system, we have 
additionally investigated the robustness of the switching 
ratio in the presence of impurity scattering. The switch- 
ing ratio acquires a higher or lower value than in the 
clean case depending, respectively, on whether the impu- 
rity configuration is blocking transport at zero magnetic 
field or not. Thus, for randomly distributed impurities 
a temperature dependent lower bound for the switching 
ratio of a sample can be set. The efficiency of magneto- 
conductance tuning then remains to be specified for the 



individual device. Conclusively, it is demonstrated that 
electron billiards of specific geometry and chains thereof 
can be used, due to regularities in the suppression of their 
transmission, to design low temperature magnetoconduc- 
tance. 
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